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ABSTRACT 

This  thesis  developed  expressions  for  the  aperature  fields  and  inte- 
gral equations  for  the  far  field  radiation  characteristics  of  the  Parallel 
Plate  Luneberg  Lens  Microwave  Antenna  operating  in  the  TEM  or  TE  _  modes. 
These  expressions  were  programmed  on  a  digital  computer  to  predict  the 
far  field  radiation  patterns  for  several  feed  systems.   Experimentation 
produced  far  field  radiation  patterns  that  were  very  close  to  the  theor- 
etical patterns  for  the  TE   mode  lens  and  substantially  different  for 
the  TEM  mode  lens. 
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I.   INTRODUCTION 

Theoretical  work  by  R.K.  Luneberg  [Ref.  1]  on  optics  in  a  medium  of 
variable  index  of  refraction  resulted  in  the  development  of  the  Luneberg 
Lens  which  has  many  applications  in  microwave  antennas.   Two  of  the  most 
common  Luneberg  Lenses  are  of  the  spherical  and  cylindrical  types,  in 
which  the  index  of  refractions  is  a  function  of  radius  only.   In  these 
lenses,  the  index  of  refi-action  is  varied  so  that  the  radiation  from  a 
point  source  placed  at  the  edge  of  the  lens  produces  a  collimated  beam 
on  the  opposite  side  of  the  lens.   The  index  of  refraction  is  1.0  at  the 
edge  of  the  lens  so  that  it  is  matched  to  free  space. 

Spherical  or  cylindrical  Luneberg  Lenses  are  constructed  by  bonding 
concentric  dielectric  shells  or  hoops  together,  or  by  a  special  mixing 
and  blending  process  which  produces  a  continuously  varying  index  of  re- 
fraction as  a  function  of  the  radius  [Ref.  2],  Cylindrical  lenses  are 
also  produced  by  taking  a  cylinder  of  dielectric  material  and  machining 
one  surface  until  the  required  average  index  of  refraction,  when  placed 
between  two  parallel  plates,  is  achieved. 

Cylindrical  Luneberg  Lenses  operate  in  either  the  TEM  or  the  TE 
modes.   If  the  thickness  of  the  lens  is  greater  than  one-half  wavelength 
but  less  than  one  wavelength,  it  operates  in  the  TE ._  mode.   However,  if 
the  lens  is  several  wavelengths  thick  it  will  operate  in  the  TEM  mode. 

Many  radar  systems  and  most  electronic  warfare  systems  require  an- 
tennas that  can  scan  a  narrow  beam  over  wide  angles  in  space.  As  long 
as  the  scan  rates  are  low,  the  present  rotating  systems  are  adequate. 
However,  when  the  required  scan  rates  are  increased  or  the  application 


requires  instantaneous  repositioning  of  the  beam,  rotating  antennas  can- 
not fulfill  the  requirements  due  to  their  inertia.   It  can  be  shown  that 
for  some  passive  electronic  warfare  systems  the  optimum  scan  rates  are  in 
excess  of  3000  rpm.   In  addition  to  this,  the  beamwidth  of  present  ro- 
tating antenna  systems  are  too  wide  to  permit  the  desired  accuracy  in  the 
determination  of  the  direction  of  arrival  of  an  intercepted  signal. 

A  new  antenna  system  is  therefore  required  to  more  fully  utilize  the 
capabilities  of  current  equipment.   The  Luneberg  Lens  has  the  potential 
to  overcome  the  above  deficiencies  and,  in  addition,  provide  a  higher 
gain  and  reduce  the  weight  of  heavy  antenna  systems  placed  on  the  super- 
structure of  ships.   In  view  of  this  potential,  this  thesis  was  designed 
to  analyze  the  parallel  plate,  continuously  varying  dielectric,  cylin- 
drical Luneberg  Lens.   The  analysis  was  to  produce  theoretical  expressions 
for  the  far  field  radiation  patterns  of  the  lens  for  any  feed  system. 
This  would  provide  a  method  of  determining  the  performance  of  any  size 
lens  with  any  feed  system  to  be  investigated  without  going  to  the  expense 
of  buying  or  constructing  a  lens  and  measuring  the  far  field  radiation 
patterns.   The  only  requirement  is  that  the  radiation  patterns  of  the 
feed  system  be  known.   The  feed  system  patterns  may  be  known  functions 
or  they  may  be  measured.   The  far  field  expressions  were  programmed  on 
the  IBM  360  computer  to  compute  and  plot  the  far  field  radiation  patterns. 
The  accuracy  of  this  analysis  was  verified  experimentally. 

Although  the  analysis  of  the  Luneberg  Lens  was  carried  out  for  arbi- 
trary feed  polarization,  lens  size,  and  frequency,  some  limits  had  to  be 
placed  upon  the  experimentation.   The  experimental  limitations  were: 
(1)  feed  systems  polarized  in  a  plane  perpendicular  to  the  plane  of  the 
lens,  (2)  the  far  field  radiation  patterns  of  the  feeds  were  used  in 


calculating  the  far  field  expressions  of  the  Luneberg  Lenses,  (3)  all  the 
experimentation  was  confined  to  x-band  due  to  availability  of  equipment, 
and  (4)  a  "roof-top"  range  was  used  for  some  measurements. 


II.   THEORETICAL  ANALYSIS 

Electromagnetic  energy  propagating  through  a  lens  follows  optical 
ray  paths  [Ref.  3].   In  the  case  of  the  radially-symmetric,  cylindrical 
Luneberg  Lens,  the  optical  ray  paths  are  such  that  the  energy  from  an 
incident  plane  wave  is  focused  into  a  cylindrical  wave  front  on  the  sur- 
face of  the  lens  diametrically  opposite  the  direction  of  incidence  as 


Figure  1.   Focusing  Properties  of  the  Luneberg  Lens. 

indicated  in  Fig.  1. 

The  refraction  necessary  to  achieve  this  focusing  effect  is  obtained 
by  properly  varying  the  index  of  refraction.   In  the  cylindrical  Luneberg 
Lens  the  index  of  refraction,  n(r)  is  given  by 

n(r)  =  /2  -(r/ro)2,  U) 
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where  r  is  any  radius  and  r   is  the  radius  of  the  lens.   Thus  the  index 

o 

of  refraction  varies  from  vl      at  the  center  of  the  lens  to  1.0  at  the 
edge. 

The  theoretical  analysis  of  the  Luneberg  Lens  is  based  upon  the 
principles  of  geometric  optics  from  which  the  Ray  Paths  and  the  Aperature 
Field  of  the  Lens  are  calculated  and  plotted.   The  far  field  radiation 
patterns  are  then  calculated  from  the  Aperature  Field. 

A.   RAY  PATHS 

The  generalized  Snell's  Law  for  circular  symmetry,  such  as  is  present 
in  the  cylindrical  Luneberg  Lens,  where  the  index  of  refraction  is  a  func- 
tion of  the  radius  only  [Ref.  3]  is: 

n(R)  sin  0  =  C   ,  (2) 

where   n(R)  =  /2-R2  =  index  of  refraction, 
R  =  normalized  radius, 

0  =  angle  between  the  ray  path  and  the  outward  meridian, 
C  =  constant. 

The  geometry  of  Fig.  2  and  equation  (2)  were  utilized  to  derive  an 
integral  formula   for  the  relationship  between  R,  Y,  and  C.  . 

(   1    4     2    2  ~1/2 
f  =  Cl  J  i   [_R  +  2R  "  Cl]    dR  '  (3) 

where  ¥  is  the  angle  between  the  direction  of  arrival  of  the  plane  wave 
and  the  radius  through  a  point  on  the  ray.   A  unique  ray  exists  for  each 
value  of  C   (0.0  <_   C   <  1.0). 


For  complete  derivation  see  Appendix  A. 
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Figure  2.   Geometry  for  Luneberg  Lens  Ray  Path  Derivation, 
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Equation  (3)  was  separated  into  two  intervals  as  given  by  equations 

(41)  and  (42)  of  Appendix  A  and  programmed  for  the  IBM  360  computer.   The 

ray  path  plots  are  given  in  Fig.  3  for  ten  values  of  C.  .   Note  that  3 

1  o 

was  equal  to  T  and  that  the  closest  point  of  approach  to  the  origin  by 
a  ray  (r  )  occured  for  0  equal  to  90°.  Each  ray  is  symmetric  about  the 
point  (r  =  r  ,  0  =  90°). 

B.   APERATURE  FIELD 

The  relationship  between  the  Feed  Distribution  and  the  Aperature 
Field  of  the  Luneberg  Lens  was  derived  utilizing  the  principles  of  geo- 
metric optics,  for  the  geometry  and  variables  defined  in  Figure  4.   The 
energy  flow  into  a  Differential  Element   dBdx  must  be  identical  to  the 
energy  flow  out  of  the  Differential  Element  dxdy  t    since  from  the  prin- 
ciples of  Geometric  Optics  it  has  been  shown  in  section  II-A  that  energy 
flow  is  along  Ray  Paths.   If  a  tube  of  rays  is  considered,  then  the  energy 
flow  is  constant  through  various  cross-sections  of  the  tube.   Therefore, 
the  energy  flow  into  the  wedge-shaped  tube,  d3dx  ,  at  the  feed  is 
S  (3,x)d3dx  ,  where  S  (3,x)  is  the  energy  distribution  per  unit  angle  in 
the  primary  feed  pattern. 

Equating  the  energy  that  enters  the  flux  tube,  d3dx  ,  to  the  energy 

leaving  the  flux  tube,  dxdy  ,  S  (x,y)  dxdy  =  S  (x,3)  d3dx  ,  where  S  (x,y) 

■'     a  J.  a 

is  the  energy  in  the  lens  aperature.   However,  S  (x,y)  is  proportional 

3. 

to     E'       ,   where      E'    is    the  magnitude   of   the   electric    field   in    the   aper- 
1    a1  'a1 

ature.      Therefore, 

S    (x,y)dxdy     -   cjE'|2dxdy     =   S-(x,8)dBdx  (4) 

where  0,     is  a  constant.   From  the  geometry  of  Figure  4 

y  =  r  sin  H'  (5) 

o 
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Figure    3.      Ray   Paths   Through   Luneberg   Lens 
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Figure  4.   Geometry  for  the  Derivation  of  the  Aperature  Fields 
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where  r   is  the  lens  radius.   r   is  defined  as  the  point  on  a  ray  nearest 
o  m 

the  origin,  and  due  to  the  symmetry  of  the  Lens,  the  ray  path  is  symmetric 
about  the  point  r  =  r   and  hence  Y  =  $   [Ref.  4].   Then  from  equation  (5) 
and  this  fact 


"J^7  =  ~Ja    =  r   cos  Y' 
dY    dp     o 


Substituting  this  into  equation  (4)  yields 
2   B1(»,3)  d3   S^x.Y) 


E' 

a1      C„    dy    C0r  cos  V 
2  2  o 


or 


|E'|   =  C   /S  (x,H0  sec  Y  ■    (6) 

3        -5     J- 

The  Aperature  field  can  therefore  be  calculated  for  any  feed  distri- 
bution for  either  the  TEM  or  the  TE   mode  Luneberg  Lens.   Two  assumptions 
were  made.   First,  that  the  feed  energy  distribution  is  a  separable  func- 
tion of  x  and  f.  Therefore 


S1(x,H')  =  S2(x)  S3(¥).  (7) 

The  second  assumption  is  that  S  (x)  and  S _(f)  can  be  calculated  or  ap- 
proximated from  the  geometry  or  by  computer  curve  fitting  of  a  measured 
feed  pattern.   For  either  type  Luneberg  Lens  the  S   (y)  can  be  equated 
to  the  feed  distribution  in  the  plane  of  the  lens  [Ref.  4],   S  (x) ,  how- 
ever, is  different  for  each  type  of  lens  and  is  a  function  of  the  lens 
geometry.   In  the  TE._  mode  lens  as  is  suggested  from  wave  guide  theory 

[S9(x)]T„   «  cos  (f),  (8) 

i  1E1Q        a 

where  a  is  the  thickness  of  the  lens.   In  the  TEM  mode  lens  the  assumption 
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was  made  that  the  S  (x)  is  constant  over  the  x-dimension  of  the  lens,  and 

[S2(x)]TEK  -  CA.  (9) 

The  Luneberg  Lens  was  designed  to  produce  a  constant  phase  in  the 
plane  of  the  lens  along  a  line  tangent  to  the  aperature  center;  therefore, 
the  aperature  fields  have  a  phasing 

kr  (1  -  cos  y).  .  (10) 

o 

Then   for   the  Luneberg  Lens 

E«    =    |E»|    eJkro(1  "    cos   *>    .  (11) 

a        'a1 

Combining  equations    (6),    (7),    and    (9)    for   the   TEM  mode  Luneberg   Lens, 

Ef    *   /S-CY)    sec  Y    .  (12) 

a  3 

However,    S    (Y)«    (E'(Y))2,    so   that 

E'(x,Y)    =   CCE    (Y)    /sec   $      ,  (13) 

a  i»   a 

where   C     is   a   constant   of  proportionality   and  E    (Y)    is    the  magnitude   of 
_)  a 

the   feed  distribution.      Rewriting  with   the  phasing  information   of  equation 
(10), 

E'(x,Y)    =   CCE    (¥)    /sec   V  ejkro(1  "    C0S   Y) .  (14) 

a  j   a 

Combining  equations    (6),    (7),    and    (8)    for  the   TE       mode   Luneberg  Lens 

the   aperature    field   is 

E'    «    Ie'I    /S_(Y)    sec   V   cos    (-) 
a        '    a1        3  a 


or 


E'(x,Y)    =    C^    cos    (-)    E    (T)    ZsTcT  ejkr°(1   "    C0S    *\  (15) 

a  o  a       a 
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C.   ELECTROMAGNETIC  FIELD  CALCULATIONS 

The  Electric  Field  of  the  Luneberg  Lens  was  calculated  utilizing  the 
aperature  field  method  and  Huygens  Principle  as  developed  in  Chapter  2 
of  Ref.  5.   The  Electric  Field  equations  are: 


and 


where 


E  = VxF —   V  x  V  x  A. 

e         wue 


H  = 


A  = 


_  V  x  A J—  V  x  V  x  F 

y        wye 


y 
4tt 


—    c 


F  = 


4ir 


f  - 

j(x;r)- 

ff  - 

K(x',V')- 


-JkR 


R 
r-jkR 


dS'  , 

dS'  , 


(16a) 


(16b) 


(17) 


(18) 


J(x'  ,H" )  -  n  x  H   , 
a  ' 


(19) 


K(x'  ,r)  -  E   x  n  t 

Si 


(20) 


dS'  =  p  dx'd^'   , 


(21) 


2v_ 
A 


(22) 


and  the  remaining  variables  are  shown  in  Fig.  5.   Only  the  Electric  Field 
will  be  considered  from  this  point  since  the  Far-Field  is  essentially  a 
plane  wave,  where  E  and  H  are  related  by  the  intrinsic  impedance  of  Free- 
Space. 

2 
The  Far-Field  equations   derived  from  equations  (16-22)  and  the  geo- 
metry of  Fig.  5  for  a  feed  polarized  in  the  x-direction  are: 


'See  Appendix  B  for  derivation. 


X 


Figure   5.      Luneberg   Lens   Geoir.etry   Used   in   the   Derivation   of   Radiation   Ex- 
pressions . 
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E0(R,*,O) 


_  jkp    e 


jkp 


4tt 


[1   -    cos    0] [cos    $] 


(23) 


-jkjR-jJ      -jkp    cos    H" 
[- r|-i HE^x'.^cos    r]    dx'    d¥\ 


E$(R,*,0) 


jkp     ft          IklRi  I         jkp    cos   H'' 
jkp    eJ  e   J     '    1[    e   J . 


4tt 


II 


R 


|K1| 


(24) 


[E      (x',r)][sin  0   sin   H"    -   cos  -T    sin   $(1  +  cos   0)]dx'    dY' 


and  E_(R,$,0)    -   0; 

K 


(25) 


where   E   (x',f  )  is  the  magnitude  of  the  Lens  aperature  field  polarized 

tlX 

in   the   x-direction   and 


2  2  2 

R       =    [R     +    (x')      +  p     +   2R(x'cos    0  cos    *  -  P    sin   ¥' 


1/7 

sin  0  sin   $  +  P    cos   H'1    cos    0)  ]         . 


(26) 


The   far-field  equations    for   a   feed  polarized   in   the  y-direction  are: 


E   (R,$,0) 


jkp  e 


jkP 


rt        -jk|R]J      -jkP   cos    V' 


4tt 


[: 


R. 


-] 


(27) 


[E      (x'  ,*')]  [cos    ¥'    sin   *    (l+cos©  )    -   sin   0  sin   H"]dx'    dY* 
ay 


E    (R,*,0) 


-  JkP  e 


jkP 


4tt 


[cos    $(cos    0  +    1) ] 


(28) 


t        -jk|Ri|      -jkp  cos    ¥' 

[- r§ — I ][E      (x',H")cos    1"]    dx'    dV 

I  R    |  ay 


and  ER(R,$,0)    -  0; 


(29) 


where        E      (x1,^')    is    the  magnitude   of   aperature   field  polarized   in   the 
ay 

y-direction  and    [R   |    is    as    given   in   equation    (26). 
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The  Far-Field  radiation  patterns  are  identical  as  can  be  seen  by 
comparing  equations  (23)  and  (24)  to  equations  (28)  and  (27)  respectively, 
Due  to  this  identity,  the  Far-Field  radiation  patterns  are  independent  of 
the  Feed  Polarization.   Therefore,  only  one  set  of  equations  need  be 
utilized  in  calculating  and  plotting  the  far-field  radiation  patterns 
for  any  feed  polarization  in  the  x-y  plane. 

Approximating  — —  by  -—■--  and  |R  |  by  |r|  -  |R'*a  | (far  field  con- 

3  \\\       |r| 

dition)   then  equations  (13)  and  (14)  become: 


E.  ■  t-ttU  ~  cos  ©][cos  $]eJ 
O   4ttK 


e    E   (x)  dx   e  J  cos  ^E   (Y)  dY 
al        J  al 


and  (30) 


E,  = 


jkp   jk(p-R) 


$   4ttR 


e    E   (x)  dx 
al 


eCER  cos  YE   (¥) 
al 


[  sin  0  sin  f  -  cos  G  sin  $  cos  V  -   sin  $  cos  j']dY  ,    (31) 

where       CEx  =  jkx  sin  0  cos  $  (32) 

and         CER  =  jkP(cos  Y  cos  0  -  sin  ¥  sin  0  sin  $  -  cos  Y) .     (33) 

Equations  (30)  and  (31)  were  programmed  for  the  IBM-360  computer  for 
calculating  and  plotting  the  normalized  far-field  radiation  patterns. 
The  computer  programs  are  given  in  Appendix  C,  and  the  theoretical  far- 
field  radiation  patterns  are  given  in  Section  III. 


3 

Calculations  show  that  for  the  worst  condition  (Y  =  90°)  at  a  distance 

of  2L,2/\  the  phasing  error  is  less  than  .06A   (21°)# 
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III.   EXPERIMENTAL  VERIFICATION 

A.   DISCUSSION 

Experimental  verification  was  undertaken  to  determine  the  validity 
of  the  equations  that  describe  the  far-field  radiation  patterns  of  the 
Parallel-Plate  Luneberg  Lenses.   In  order  to  accomplish  this,  two  Lune- 
berg  Lenses  were  acquired.   One  was  a  TEM  mode  lens,  44  inches  (37.2  A 
at  10  GHz)  in  diameter  and  2.38  inches  thick.   The  other  was  a  TE 
mode  lens,  10.25  inches  (8.67A  at  10  GHz)  in  diameter  and  0.5  inches 

thick.   All  measurements  were  made  with  the  lenses  operating  in  the 

2 
receive  mode  and  for  far-field  conditions  (distance  >  2L  /A).   Measure- 
ments on  the  TE   mode  lens  were  conducted  in  an  anechoic  chamber  and 
those  for  the  TEM  mode  lens  were  conducted  in  a  "roof-top"  range.   The 
transmitting  set-up  is  shown  in  Fig.  6,  and  was  operated  in  the  CW  mode 
at  10  GHz.   The  only  exception  was  for  a  frequency  characteristics  check 
of  the  TE  _  mode  lens.   The  receiving  set-up  is  given  in  Fig.  7. 

The  procedure  followed  was  to  measure  the  far-field  patterns  of  the 
feed  and  the  Luneberg  Lens  and  then  using  the  feed  data,  calculate  the 
far-field  patterns  and  compare  the  two  results. 

It  was  noted  that  in  the  case  of  all  feeds  with  both  lenses  that  when 
the  feed  was  placed  very  close  to  the  lens  the  first  side  lobes  almost 
merged  into  the  main  beam.   As  the  distance  between  the  feed  and  lens 
increased  the  side  lobes  became  distinct  and  dropped  very  rapidly.   Then 
the  particular  placement  of  the  feed  was  found  which  resulted  in  optimum 
(lowest)  side  lobe  levels.   When  the  separation  was  increased  to  a  dis- 
tance greater  than  the  optimum,  the  main  beam  spread  out  and  engulfed 
the  first  side  lobes. 
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B.   TE  Q  MODE  LUNEBERG  LENS 

The  TE   Mode  Luneberg  Lens  was  fed  by  two  different  feed  systems. 
The  first  was  an  "open-end"  waveguide  feed  (hereafter  referred  to  as 
the  waveguide  feed)  and  the  other  was  a  two  element  Yagi  array  (driven 
element  plus  director)  resonant  at  10  GHz.   The  far  field  patterns  of 
the  lens  aperature  fields  for  these  feeds  are  given  in  Figures  8  and  16 
respectively.   Fig.  9  is  a  comparison  of  the  theoretical  vs.  the  experi- 
mental far  field  patterns  in  the  planes  of  the  lens  (E-^plane),  when  fed 
by  the  waveguide  feed  at  a  frequency  of  10  GHz.   The  gain  over  isotropic 
of  this  system  is  approximately  20  db,  as  compared  to  a  16  db  standard 
horn,  with  a  one-half  power  beam  width  of  6.5  degrees.   The  salient 
points  of  the  two  curves  are  given  in  Table  I,  below. 

TABLE  I 

3  db      First  Null   First  Sidelobe   First  Sidelobe 
Beam  Width   Location     Location         Amplitude  (db) 

left  right  left  right  left  right 

Theoretical     6.8°     -8.2°  8.3°  -10.8°  10.6°  -19.6  -19.6 

Experimental    6.5°     -7.4°  7.8°  -  9.9°  10.3°  -17.9  -18.2 

Difference      0.3°      0.6°  0.5°  0.9°   0.3°  1.7    1.4 

Figure  10  is  a  comparison  of  the  theoretical  and  experimental  H-plane 
(perpendicular  to  plane  of  lens  at  beam  center)  far  field  patterns  for 
the  waveguide  feed. 

In  order  to  determine  the  frequency  characteristics  of  the  TE   mode 
lens  when  fed  by  the  waveguide, feed  far  field  patterns  for  7,8,9  and  11 
GHz  were  also  calculated  and  are  presented  in  Figures  11  and  15.   The 
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experimental  far  field  patterns  were  different  from  the  theoretical  pat- 
terns by  about  the  same  amount  as  the  patterns  at  10  GHz. 

The  experimental  vs.  theoretical  E-plane  far  field  patterns  of  the 
TE   Mode  Lens  fed  by  the  Yagi  feed  are  given  in  Fig.  17.   The  differ- 
ences in  side  lobe  and  null  locations  and  side  lobe  levels  are  negligible, 
The  gain  over  isotropic  of  this  arrangement  is  about  21.5  db  as  compared 
to  the  16  db  standard  gain  horn,  with  a  one-half  power  beam  width  of 
6.2  degrees. 

A  comparison  of  the  lens  aperature  fields  (  Figures  8  and  16)  show 
that  the  Yagi  array  feed  utilizes  the  lens  aperature  more  efficiently 
than  the  waveguide  feed. 

C.   TEM  MODE  LUNEBERG  LENS 

The  theoretical  and  experimental  E-plane  far  field  patterns  of  the 
TEM  Mode  Luneberg  Lens  when  fed  with  a  dielectrically  loaded  sectoral 
horn  feed  (feed  pattern  and  lens  aperature  field  given  in  Fig.  19)  are 
given  in  Figures  20  and  21,  respectively.   The  difference  between  these 
two  patterns  is  about  1.0  degrees  in  the  one-half  power  beam  width,  0.5 
degrees  in  the  first  null  placement,  2.5  degrees  in  the  first  sidelobe 
placement,  and  8.0  db  in  first  sidelobe  amplitudes.   The  gain  over  iso- 
tropic of  this  arrangement  was  about  25  db .   The  theoretical  H-plane 
far  field  pattern  is  given  in  Fig.  22.   The  3  db  theoretical  H-plane 
beam  width  is  22  degrees  as  compared  to  25.5  degrees  obtained  experi- 
mentally by  the  Naval  Electronic  Laboratory  Center  [Ref.  2]. 

When  the  TEM  mode  lens  was  fed  by  the  waveguide  feed  (Fig.  8)  the 
theoretical  (Fig.  18)  and  the  experimental  far  field  pattern  differences 
were  somewhat  less  that  they  were  for  the  sectoral  horn,  but  due  to 
range  difficulties    reliable  experimental  patterns  were  not  obtained. 
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It  is  felt  that  the  differences  between  the  theoretical  and  experi- 
mental patterns  were  too  great  to  warrant  further  investigation  without 
refinement  of  the  feed  systems  analysis. 
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Figure  8.   Waveguide  Feed  Pattern  and  Luneberg  Lens  Aperature  Field 
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Figure  9.   Luneberg  Lens  Far  Field  Radiation  Patterns  with  Waveguide  Feed, 
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Figure  10.   Luneberg  Lens  (TE   mode)  H-plane  Far  Field  Patterns  for 
Waveguide  Feed. 
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Figure  11.   Luneberg  Lens  E-Plane  Patterns  for  Waveguide  Feed  at  7.0  GHz 
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Figure  12.   Luneberg  Lens  E-Plane  Patterns  for  Waveguide  Feed  at  8.0  GHz 
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Figure  13.   Luneberg  Lens  E-Plane  Patterns  for  Waveguide  Feed  at  9.0  GHz 
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Figure    14.      Luneberg   Lens    E-Plane  Patterns    for  Waveguide   Feed   at    10   GHz, 
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Figure    15.      Luneberg  Lens    E-Plane   Patterns    for  Waveguide   Feed   at    11   GHz 
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Figure  16.   Yagi  Array  Feed  Pattern  and  Luneberg  Lens  Aperature  Field  for 

Yagi  Feed. 
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Figure   17.      Luneberg   Lens    (TE        mode)    E-Plane   Patterns   For  Yagi   Feed   at 
10   GHz. 
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Figure  18.   TEM  Mode  Luneberg  Lens  Theoretical  E-Plane  Radiation  Patterns 
With  Waveguide  Feed  at  10  GHz. 
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Figure  19.   Dielectric  Loaded  Sectoral  Horn  Feed  Patterns  and  Aperature 
Field  for  Sectoral  Horn  Feed. 
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Figure  20.   TEM  Mode  Luneberg  Lens  Theoretical  E-Plane  Radiation  Patterns 
For  Sectoral  Horn  Feed  at  10  GHz. 
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Figure  21.   TEM  Mode  Luneberg  Lens  Experimental  E-Plane  Pattern  for 
Sectoral  Horn  Feed  at  10  GHz. 
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Figure  22.   TEM  Mode  Luneberg  Lens  H-Plane  Radiation  Patterns  for  Sectoral 
Horn  Feed  at  10  GHz. 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  the  experimentation   did  not  show  conclusively  the 
approach  of  this  thesis  to  the  problem  of  predicting  the  far  field  radia- 
tion patterns  of  the  Luneberg  Lens  is  valid.   However,  the  results  did 
indicate  that  the  method  is  valid  for  the  TE  n  mode  lens.   There  are  at 
least  two  reasons  for  the  differences,  assuming  that  the  approach  is 
valid. 

First,  the  far  field  patterns  of  the  feed  systems  were  used  in  cal- 
culating the  aperature  fields  and  far-field  patterns  of  the  lens,  when 
in  fact,  the  lens  was  not  in  the  far  field  of  the  feeds  except  for  the 
case  of  the  Yagi  array  on  the  TE   _  lens.   This  resulted  in  errors  in 
the  aperature  fields  used  in  calculating  the  far  field  patterns. 

The  results  of  the  TE   mode  lens  seems  to  confirm  that  this  was  in 
fact  one  of  the  major  sources  of  error.   The  far  field  of  the  waveguide 
feed  is  approximately  4  cm.  and  it  was  operated  at  1.7  cm  from  the  lens 
for  optimum  side  lobe  levels.   The  differences  between  the  theoretical 
and  experimental  patterns  were  an  average  of  0.6  degrees  in  first  side- 
lobe  placement  and  1.5  db  in  first  sidelobe  amplitude.   However,  when  the 
Yagi  array  feed,  with  a  far  field  distance  of  0.4  cm.,  was  used  at  a  dis- 
tance of  1.7  cm  from  the  lens,  the  differences  in  first  sidelobe  place- 
ment and  amplitude  were  negligible. 

The  far  field  of  the  dielectric  loaded  sectoral  horn  is  approximately 
9.6  cm.   It  was  operated  at  about  2  cm  from  the  TEM  mode  lens  and  the 
results  were  different  by  about  8  db  in  the  first  sidelobe  levels  and 
2.0  degrees  in  the  first  sidelobe  placement.   When  the  TEM  mode  lens 
was  fed  with  the  waveguide  feed  the  first  sidelobe  differences  were  on 
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the  order  of  4.0  db  in  amplitude  and  1.0  degrees  in  placement.   The  far 
field  patterns  of  the  TEM  mode  lens  with  the  waveguide  feed  were  not  pre- 
sented in  section  III  due  to  the  lack  of  confidence  in  the  results. 

Second ,  the  assumption  was  made  that  the  energy  distribution  over  the 
x-dimension  of  the  TEM  mode  lens  is  constant,  while  in  fact  it  cannot  be 
due  to  the  metal  plates  at  the  top  and  bottom  of  the  lens.   An  attempt 
to  measure  the  aperature  field  distribution  over  the  x-dimension  of  the 
lens  failed  due  to  the  disturbance  of  the  field  caused  by  the  dipole 
probe  used.   The  x-variation  of  the  aperature  field  did,  however,  appear 
to  be  a  standing  wave  between  the  two  metal  plates. 

It  is  the  opinion  of  the  author  that  the  approach  of  this  thesis  to 
the  far  field  analysis  of  the  Luneberg  Lens  is  valid,  and  that  the  dif- 
ferences between  the  theoretical  and  experimental  results  is  almost  com- 
pletely due  to  the  feed  pattern  information  used  in  calculating  the 
theoretical  patterns. 

In  line  with  the  above  thoughts  the  following  recommendations  for 
future  investigations  are  made: 

A.  That  the  near  field  (theoretical  or  experimental)  of  the  feeds 
be  used  in  the  calculation  of  the  far  field  radiation  patterns  of  the 
Luneberg  Lenses. 

B.  That  a  broad-band  Yagi  array,  such  as  Yagi  constructed  from 
capacitative  lumped  loaded  cylindrical  elements,  be  used  to  feed  the 
Luneberg  Lens. 

C. •  That  multiple  Yagi  arrays  be  used  to  determine  the  possibility 
of  using  electronically  scanned  Yagi  array  (such  as  in  item  B)  to  feed 
the  Luneberg  Lens  in  a  direction  finding  system. 
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D.   That  the  computer  programs  of  Appendix  C  be  modified  to  calculate 
the  far  field  patterns  of  the  Constant-K  Lens.   This  would  require  a 
change  of  function  subprograms  CFCT  3  and  CFCT  4  to  reflect  the  geomet- 
ric optics  of  the  Constant-K  Lens. 


APPENDIX  A 
RAY  PATH  CALCULATIONS 

A.   DERIVATION  OF  RAY  PATH  FORMULA 

The  generalized  Snell's  Law  for  circular  symmetry  is  given  in  Ref.  3 
as 

n(R)R  sin  0  =  C  (34) 

when 

n(R)  =  Index  of  refraction 
R  =  Normalized  radius 

0  -  Angle  between  ray  path  and  outward  meridian 
C  =  Constant 

From  the  geometry  of  Figure  23, 


sin  0  =  A*"1        =       R  (35) 

./mr>2  +  (dR)2       /r2  +  (§) 2 

Substituting  equation    (35)    into    (34)   yields 

'    nR2  CldR 

Ci=  or      d¥  =  X  .  (36) 

A     +   (^)  R/n   R     -   C2 

For  the   Luneberg  Lens    (normalized   radius), 


n  =    fa   -   R2  (37) 

which  yields 

C   dR 
d¥  = (38) 


r/rS- 


or 


2R     -    C 


I  |[-R4  +  2R2  -  Cx2]  1/2dR.  (39) 


46 


Figure  23.   Geometry  For  Luneberg  Lens  Ray  Path  Derivations 
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Integrating  this   equation 

1  -1     /        2        1  -12        2  2  /        2 

f  -     j  sin        /l-C1Z  -  |  sin      [(R-Cp/(RVl-Cp].  (40) 

This  equation  gives  four  values  of  4*  for  each  value  of  C  and  R  and 

hence  two  rays.   Geometrical  considerations  produce  the  following  for- 
mulas: 

Vx  -  j  -  sin"  A-C   "  +  V  (41) 


and 


where 


*2  =  7T  -  Y  ,  (42) 

R      c 

V  =  f dR  , 


1  R  /-R4  +  2R2  •  C  2 


4'   corresponds  to  4*  for     4*   <  4'  <  tt  , 
1  m  —   — 


4>   corresponds  to  Y  for     4'   <  4'  <  4' 
2  o        m, 


and  4/  is  the  value  of  4'  for  which  R  =  r  . 

m  m 


B.   CALCULATIONS  OF  r  . 

m 

For  each   ray  path    the   closest   point   of  approach    to   the   origin,    r    , 
occurs    for   0  =   90°.      Therefore,    from  equations    (34)    and    (37) 

r    [n(r   )]    sin   90°   =   C, , 

mm  1 

or 

r   A-x   2      =      C,. 

mm  1 

Squaring  both  sides  yields 

4      2     2 
r   -  2r   +  C,   =  0 

m      m     1 
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or,  by  the  quadratic  formula, 


r  -  1  ±  /l-c,2  . 

m  1 

The  negative  square  root  must  be  taken  since  the  positive  square  root 

gives  an  r   greater  than  the  normalized  radius.   Therefore, 


m 


r  -  (1  -  £7*)  1/2 


(1  -  /1-C  )  x/"  ,  (43) 

m  1 

which  is  the  closest  point  of  approach  of  the  ray  path  to  the  center  of 

the  lens . 
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APPENDIX  B 
ELECTROMAGNETIC  FIELD  DERIVATION 

A.   DERIVATION  OF  ELECTRIC  FIELD  EQUATIONS. 

The  field  equations  for  any  aperature  on  which  the  electric  and  mag- 
netic currents  are  known  can  be  calculated  utilizing  the  Aperature  Field 
method  and  Huygen's  Principle  as  developed  in  chapter  2  of  Ref.  5.   The 
electromagnetic  field  equations  for  any  aperature  [Ref.  5],  and  the  geo- 
metry of  Figure  24  are: 


and 


where 


E  = yx  F —  V  x  v  x  A 

e         wye 


11=   -  V  x  H 

y         wye 


V  x  V  x  F; 


4rr 


ff  J(x',Y')e-jklRll 


R, 


F  = 


4tt 


f  K(x',V)e  jk'Rl 


J(xf,T')  =  n  x  H 


dS 


dS' 


(44) 

(45) 
(46) 

(47) 
(48) 


K(xf ,Y')  =  E  x  n 


(49) 


^ 


(50) 


dS'  =  P  dx'  dY' 


(51) 


2      2 


1/2 


1^1=  [R+(R')  -2R(x,cos$cos0-psiny'sinesin$+pcosH"cos0)]  ; 

(52) 

n  =  -sin  4" a  +  cos  V'a  =  unit  normal  to  aperature.     (53) 
y         z  ^ 
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$   a 


Figure  24.   Luneberg  Lens  Geometry  Used  in  the  Derivation  of  Radiation 
Expressions . 
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Equation  (1)  was  divided  into  two  parts: 


where 


and 


E  "  El  +  E2 


E   = v  x  F 

1      E 


E   = d-   V  x  V  x  A, 

2     wue 


and  considered  separately. 

From  equations  (47),  (49)  and  (54) 


(53) 


(54) 


(55) 


1  e  4tt 


k(x'.t')  e-3k|Rilds.] 

Ir,  I 


(56) 


Taking   all   constants   outside   the   curl  equations    and   interchanging   the 
order   of   the   integration  and   curl   operators, 


El  4tt 


V   x    [Mx    >]    }? ]    dS\ 

R-. 


(57) 


Using   the  vector  identity 


V   x   Kf  =  Vf   x   K  +   f  (V   x   K) 


where 


rJklRil 


f  = 


R, 


equation    (57)    reduces    to 


(58) 


(59) 


El        4tt 


since 


1       e"JklRli      -        - 
(jk  +TC"]-)(-    |R   |       »ar  ><   K(x',f'))    dS'    , 


vf     ptt      JL_)e-dHRil 


\\       \\\ 


and 


V    x   K  =   0, 


K  is   independent    of   the  point   of   observation. 
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(60) 


From  equations  (46)  and  (55) 


E0  =  -J—  y  x  V  x  [y^ 


J(x',Y ')e  J  '  n 


ds'] 


R, 


Taking  all  constants  outside  the  curl  equations  and  interchanging  the 
order  of  integration  and  the  curl  operators 


_  "J 


2   4 


iTwe 


r  Tt    t   i  \  -jkl  Ri  I 

—   —    rJ(x  ,y  )e  J  '  -1-1 


V  x  V  x  [ 


R, 


]  ds'. 


(61) 


Applying  the  vector  identity  of  equation  (58)  to  the  integrand  twice 
yields 

2 


V  x  V 


,  „  (J(*'>^-3k|Rli)  .  (-  JL_  +  Jttj  +  _2_)e-JklRll  [a_,(I_xJ) 


R. 


1*1 1    lR! 


R. 


r 
(62) 


Then   from  equations    (61)    and    (62) 


E     =        J 
2        47Tooe 


|R, 


2          e"jklRll  -         - 

-^-y)(- )k    x   (a      x  J)]    ds'. 

|R,  I  IRJ 


(6  3) 


Equations    (51),    (53),    (60)    and    (63)    yield   the   Electric   Field      or   any 

cylindrical   aperature    field   at   any   distance    from   the   aperature.      For   the 

2  i        3 

far   field  only,    the   1/  |R    |      and   1/  |R    |      terms    can  be   neglected   leaving 


E      =  ^ 

1  4  TV 


rre-jklRil  -      _ 


(a      x  K(x'  ,  V'))dx'    d¥' 


and 


2       ,,    -jk  |R     | 

—[a    x     (a    xJ(x',f)]dx'    dV  . 

I*!  I 


E9     =     i 

£        h  ttooe; 


JJ 


(64) 


(65) 


The  Luneberg   Lens  was    designed   to  produce   constant   phase   along   a    line 
perpendicular   to   the   aperature   center  and   therefore   the   aperature    field 
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has   a  phase   of  kp(l  -   cosy1)      so   that: 


E     =   E        e 
a  a1 


jkp(l  -    cos   4") 


and 


H     =   H        e 
a         a. 


jkp(l  -   cos   y1) 


E        and  H        are   the   E  and  H   fields   of   the   aperature 
al  al 

K  and  J    can   then  be   rewritten   as 


and 


K  -    (E     x  n)ejkp(1  "   COS   *'>   =  Fe^P(1  "   COS   *'> 

al  X 


T        /-     5-    \    JkP(l   -    cos    4")        -     jkP(l  -    cos    ¥') 
J  =    (n   xR      )eJ  =  J.eJ 

a,  i 


(66) 


(67) 


Equations    (64)    and    (65)    can  be   rewritten   as 


and 


-         jkP    ((  e-^klRllejkP(1-   COS   r)  = 


El    ~    ATI 


E     =  Jk   p 

2        Attojc 


R. 


F      dx'    d*i" 


-jk|Ri|    jkP(l  -   cos   y') 
(- — )F      dx'    d«P' 

IrJ 


(68) 


(69) 


where 


and 


F.  =  a  x  K(x'  il>')  =  a   x  (E   x  n)  , 
1    r       ,r      r     a 


(70) 


F  =  a  x  (a  x  J  (x',y')  =  a   x  [a   x  (n  x  H   )]      (71) 
From  equation  (70)  and  vector  identity 


A  x  (B  x  C)  =  (A  •  C)B  -  (A  •  B)C, 


(72) 


or 


where 


F  =  (a  •  n)E   -  (a   •  E   )n. 
1     r      a      r    ai 


Fl  =  C7Ea  -    V  ' 


C  -  a   •  n, 
7    r 


(73) 
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ana 


C_   a   •  E   . 
8    r    a 


(75) 


Assuming  that  the  aperature  field  is  polarized  in  the  x-direction  (i.e, 


E  -  E   a   and  H   =  H   a  ) 


ai 


a±  ay  y 


C  =  -  sin  0sin  $sin  ¥'  +  cos  Ocos  Y' 


(76) 


and 


Cn  =  E  sin  Gcos 
8    ax 


(77) 


Combining  equations  (73),  (76)  and  (77)  gives 


F,  =  (-sin  0  sin  $  sin  Yf  +  cos  0  cos  Y)E   a 
1  ax  x 

+  E   sin  0  cos  $  sin  ¥'  a 
ax  y 

-  E   sin  0  cos  $  cos  Y1  a 
ax  z 

Converting  to  spherical  coordinates 


(78) 


where 


F,  =  F,   a  +  F„rt  art  +  F  .  a 

1    lr   r    10  0    i$   $ 


F.   =  0 
lr 


F,„  =  E   cos  T'  cos  $  , 
10    ax 


(79) 

(80) 
(81) 


and 


F„ ,  =  E   (sin  0  sin  7'  -  cos  0  cos  Y1  sin  $) 
1$    ax 


(82) 


When   the   above   equations   are   substituted  into  equation    (68) 


e  -Jte. 

4tt 


/"  e-^^lle^^1  "   cos    H") 


[F-    a     +  F.a     +  F.    a    ]dx'    dV    (83) 

lr  r  10  0         1$  $ 


or 


E.   +  E.    a     +  E_    a     +  E.    a 
1-         lr   r  10  0  1$   $ 


(84) 


where 


E,      =  0, 
lr 


(85) 
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1$        4tt 


/   p-JklRll    J1^1  ~    cos   ^') 
F,,    dx'    d4" 

R  J 


(86) 


and 


f        =  Jkp_ 
10        4tt 


/■      -jk|  R]J    jkp(l  -   cos   4'') 


R. 


F        dx'    dY' 


(87) 


Similarly,    from  equations    (71)    and    (72) 


F     =    a      x    (CQn   -    C_H   ) 
2  r  9  7  a 


or 


where 


F2   =   S(ar  X   n)    "   C7Ur  X   V ' 


CQ   5   a      •    H 
9  r  a 


(88) 


(89) 


and 


a     x  n  =    (sin  0   sin   $   cos   ¥'   +  cos   0   sin  f')a 
r  x 

-  sin  0   cos    $   cos   ¥'    a 


sin   0   cos   $   sin  4*'    a 

z 


(90) 


Again  assuming  an   aperature    field  polarized   in   the  x-direction 


H        =  H        a 
a±  ay     y 


a     x   H     =   -  H        cos   0a     +  H        sin  0   cos    $   a    , 
r  a  ay  x  ay  z 


(91) 
(92) 


Crt   =   H       sin  0  sin   $    ,  (9  3) 

9  ay 

and   C     is   as   given  by   equation    (76).      From  equations    (76)    and    (88-93) 

—  2  2  2  — 

F_   =  H      (sin  0  sin   $   cos   ?'    +  cos   0    cos   4") a 
2  ay  x 

2  — 

-  H     sin  0  sin   $   cos    $   cos    4"    a 

ay  y 

-  H     sin  0   cos    $   cos   0   cos    4,t    a    .  (94) 

ay  z 

Converting   to  spherical   coordinates 


F_    =   F0    a     +  F_.a_   +   F. . a . 
2  2r   r  20  0  2*   * 


(95) 
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or 


ik2D    fl   e-JklRllejkp(1  "    C°S   f,)  - 

E.  =  ^-P-  (F2rar  +   F2Qa0  +  F^dx'    d¥  '  (96) 


2        4ttu)£   J 
vhere 


R. 


and 


F2r  =  0, 


F_   =  H        cos   0   cos   Y'cos 
29        ay 


F_   =  -H        sin  $    cos   Y ' 

2$  ay 


(97) 

(98) 

(99) 


When   the   above   equations   are   substituted  into  equation    (69) 


E     =  jk  P 

2        4ttu)£ 

or 


where 


-jklRjJ    jkp(l  -   cos   V) 


■][F2rar+F20ae  +  F2$a$]dx,dr 


and 


E_   =   E.    a     +  E_    a     +  E0_a    • 
2  2r   r  20  0  2$   $' 


E-      =   0, 
2r  ' 

..2 

E        =  Jk  p 

20  47TOJC 


.-2 

E        =   Jk   P 

2$        4ttco£ 


-jklR.  I     ikp(l  -   cos   Tf  ) 
j        '    1' eJ 


^-jkJR-j^^jkpd  -   cos   Y') 


(100) 

(101) 

(102) 

F20   dx'dl1"    i 

(10  3) 

FOA   dx'dY' . 
2$ 

(104) 

From  equation  (53) 


where 


E  =  E,  +  E0=  E  a  +  E  a  +  E  a 
1    2    r  r    0  0    $  < 


E  =  E.   +  E. 
r    lr    2r 


and 


E0  =  E10  +  E20> 


E.  =  E..  +  E__. 
$    1$    2$ 


(105) 

(106) 
(107) 
(108) 


Substituting  equations  (85)  and  (102)  into  (106),  (87)  and  (103)  into 
(107),  and  (86)  and  (104)  into  (108), 


E  =  0, 
r 
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E      =^ 
0        4tt 


f     e-JklRlUMl  -   cos   «p») 


■]  [E        cos    y'    cos    $  + 


ax 


—  H        cos   0   cos   y'cos    $    |    dx1    dy'    , 
toe      ay 


and 


;       =    3kP 


-jklRjl    jkp(l  -   cos   y') 


[E      (sin   0  sin   y'    - 
"    ax 


-   cos    o  cos    y'sin   $)   H H        sin    $   cos    y1 ]    dx'    dy' . 

toe   ay 

Considering  the  far  field  relationship  of  E  and  H  at  the  aperature  of 
the  lens 

H    -/I  I 

ay    y    ax 

which  when  substituted  in  the  above  equations   yields 


E     =  0, 
r 


(109) 


and 


E      =  J'--p-    [cos    $    (1  +   cos    0]e   jkp 
0        4tt 


-jk|RjJ    jkp    cos   y' 


E        cos   y'    dx'    dy' 
ax 


•  i          -i        r  r        -ik  I  R_  I    jkp    cos    y' 
E      -JkP    ejkp  [e         !    lle ][E      ] 


(110) 


[sin  0   sin   y'    -   cos    y'sin   $    (1  +  cos   o)  ]    dx'    dy'.       (Ill) 

A  similar   derivation   for  an   aperature   field  polarized   in   the  y-dir- 
ection  yields 


E     =  0, 
r 

.,         .,  n    tt     -jklR,  I    -jkp    cos   <F' 

En  -  f P-  e^kP  llg [E      ] 

0        4tt  J  J  ,R    ,  ay 


(112) 


[cos   y'sin   $    (1  +  cos   0)    -  sin  0  sin  y1]    dx'    dy',       (113) 


and 


E„   "  f^Uos    $(cos    0  +   D] 


[E        cos   V\    dxf    d«F', 
ay  58 


f     -jklR,  I    -jko    cos   vj»« 
f    c   J     '    I1 e 


l*!l 


(114) 


where  E   is  the  magnitude  of  the  aperature  field  polarized  in  the  y- 
ay 

direction. 

The  two  sets  of  equations  (109-111)  and  (112-114)  when  combined  will 
produce  the  Far-Field  expressions  for  any  aperature  field  polarized  in 
the  x-y  plane.   If  one  is  interested  in  the  Far-Field  patterns  only, 
then  comparison  of  equations  (110)  and  (111)  to  equations  (113)  and  (114), 
respectively,  indicates  that  the  normalized  Far-Field  patterns  are  iden- 
tical for  any  aperature  field  polarized  in  the  x-y  plane.   Therefore, 
in  order  to  obtain  the  Far-Field  patterns,  only  one  set  of  the  equations 
need  be  considered. 
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APPENDIX    C 
LUNEBERG    LENS    RAY    PATHS 

IMPLICIT    REAL*3( A-H,R,P) 
C  VARIABLES    USED    IN    THE     PROGRAM    ARE: 

C  CL-RAY    PATH    CONSTANT. 

C  RMIN-CLOSEST     POINT    OF    APPROACH    OF    RAV    PATHS    TO 

C  LENS    CENTER. 

C  YRAY-    Y-COMPONENT    OF     A    RAY    PATH    POINT. 

C  XRAY-    X-COMPONENT    OF     A    RAY     PATH    POINT. 

C  R-    NORMALIZED    LENS    RADIUS. 

REAL    LABEL     /•  «/ 

REALMS  TITLEU2) 
C      READ  TITLE  TO  BE  TYPED  AT  BOTTOM  CF  GRAPH 

READ(5, 100)  TITLE 
LOO    FORMAT (6A8) 

C      FCT5  IS  A  FUNCTION  MICH  WHEN  INTEGRATED  DESCRIBES 
C      THE  RAY  PATHS  THROUGH  THE  NORMALIZED  LUNEBERG  LENS. 

EXTERNAL  FCT5 

DIMENSION  XRAYC 1100), YRAYC  i 100) , AUX(25) 

COMMON  CL 
C      INITILIZE  CI  FOR  THE  1ST  RAY. 

C1=0.9D0 

PI=3. 1415926500 
C      ENTER  LOOP  TO  OBTAIN  DIFFERENT  RAYS 

DC  1  1=1,9 
C      CALCULATE  RMIN 

RMIN=DSQRT( 1 .ODO-DSQRT( 1.0D0-C1**2) ) 

R=1.0D0 

ASTP  =  0.002  5D0 

NUM=( (R-RMIN)/ASTP)+1 

RHC=DARSIN(C1  ) 

N=2*NUM 

K=N 

M  =  l 

PH=PI/2.0D0-DAPSIN(DSCRT( 1.0D0-C1**2) ) 

R  R  =  R 

X1=R*DCCS(RH0) 

3  IF( 1.0. LE. XI)     GO    TO    2 

C  CALCULATION    OF    RAY    EXTENDING    IN    FRONT    OF     LENS. 

K  =  K+1 

Xl=Xl+0.005 

XRAY(K)=X1 

YRAY(K) =C1 

GO    TO    3 
C  CALL     INTEGRATING    SUBROUTINE    TO    SOLVE     INTEGRAL    EQ 

2  CALL    DQATR(R^T1 .00,0.0000001,15 ,FCT5, DAL,  IER7, AUX) 

CALCULATE    X    AND    Y    VALUES    TO    BE     USED    IN    PLOTTING    RAYS 

PFI=PH+DAL 

PH2=PI-DAL 

XRAYLN  )=RR*DCGS(  PHI  ) 

YRAY(N) =RR*DS  IN( PHI  ) 

XRAY( M) =RR*DCOS< PH2) 

YRAY(M)=PR*DSIN( PH2  ) 

IF(N-M. LE. 1)     GO    TO    4 

N  =  N-1 

M  =  M  +  1 

RR=RR-ASTP 

GO    TO    2 
C  ENTER     INTO    CONTROL     FOR    PLOT  IMG    RAYS 

4  IF( I.GT.L )G0    TO    5 

6  M0D=1 
GC    TO    7 

5  M0D=2' 

7  CALL  ORAW(  K,  XRAY, YRAY, MOD, 0, LABEL, TITLE, 0,0,4, 4,2  ,  2,8 
18  ,0,LA5T) 

IF( MOD.NE.  1)G0  TO  300 
M0D=2 
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GO  TO  400 
300    IF( I.GE.9)  GO  TO  350 

GC  TO  400 
350    MCD=3 
400    DO  500J=1,K 
500    YRAY( J)=-YRAY(J ) 

CALL  DRAW(K,XRAY,YRAY,MOD,0  .LABEL, TITLE, 0,0, 4, 4, 2, 2, 8, 
18,0, LAST) 

Cl=Cl-0.1D0 
I      CONTINUE 

STOP 

END 

FUNCTION  FCT5(RAN) 
C      FCT5  IS  A  CUNCTION  WHICH  WHEN  INTEGRATED  DESCRIBES 
C      THE  RAY  PATHS  THROUGH  THE  LUNEBERG  LENS. 

IMPLICIT  REAL  ~8(A-H,R,P) 

COMMON  CI 

FCT5=Cl/(RAN*DSQRT(-( RAN* *4)+2.0D0*(RAN**2)-(C 1**2)  )  ) 

RETURN 

END 
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FEED  PATTERN  FITTING 
AND 
APERATURE  FIELD  CALCULATIONS 

C      THIS  PROGRAM  WAS  DESIGNED  TO  FIT  A  LEAST  SCLARE 

C      PCLYNOMINISL  TO  A  NUMBER  OF  POINTS  FOR  A  MEASURED 

C      FEED  PATTERN  AND  THEN  PLOT  THE  FEED  FUNCTION  BY  USING 

C      TFREE  TIMES  AS  MANY  POINTS  (K2)  IN  ORDER  TC  DETECT  ANY 

C      EXCURSION  THAT  THE  PCLYNOMINISL  FIT  MAY  HAVE  TAKEN 

C      BETWEEN  TWO  ADJACIENT  PCINTS  . 

IMPLICIT  REAL*8(D) 

DIMENSION  DTHETA(100),DMAG(100),DA(20),DWI(80),DY(80), 
IDELY(80),A(20) ,DSA(20)  ,THETA(200)  ,EMAG(20C)  ,DTHE(  LOO) 
2,ENA(200) 

REAL-8  DTITL(  10)/10*»  '/ 

RE/SL    LABEL/'  •/ 

REAL*8  TITLE(12) 
C      READ  GRAPH  TITLE 

PEAD(5T  100)  TITLE 
100    F0RMAT(6A8) 

C      M  IS  THE  DEGREE  OF  PCLYNOMINIAL  TC  BE  FITTEC 
C      N  IS  THE  NUMBER  OF  POINTS  IN  THE  CURVE  TO  BE  FITTED. 

READ(5,20C)M,M 
200    FORMAT( 13/12) 

C      READ  IN  DATA  POI NTS-DT HET A  IS  THE  ABSISSIA  VALUE  AND 
C      DMAG  IS  THE  SQUARE  OF  THE  ORDINATE  VALUE-SCUARE  LAW 
C      DETECTOR  WAS  USED  TO  OBTAIN  DATA. 

READ(5,400)  (DTHETAd  )  ,DMAG(  I)  ,  I  =  1,N) 

DC  60  J=1,N 
60     DNAG( J)=DSQRT(DMAG(J)  ) 
400    FCPMAT(F10.7, F12  .6) 

DO  50  1=1, N 
50     DTFE( I )=DTHETA( I)*(2.14159  2  65DO/180.DO) 

DO  5  1=1, N 

5  CM(I)  =  l.ODO 

C      LSCPL2  WILL  CALCULATE  THE  LEASE  SQUARE  PCLUNOMINIAL 
C      CCEFFCIENTS. 

CALL  LSQPL2(N,M,DTHE   ,  DM  AG ,DW I , DY , DEL Y , DA  ,  D SA , DT  ITL ) 

L=1+IABS(M) 

DO  6  1=1, L 

6  A(  I  )  =  DA( I ) 
WRITE(6,300)  (I, A(  I)  ,1  =  1, L) 

300  FOPMAT(«  •  ,  «A( •  , 14, « ) =«  ,F18.6) 
THET=DTHE( 1) 

K2=3*N 
C      CALCULATE  SIZE  OF  STEPS  TO  BE  USED  IN  PLOTTING  CURVE. 

STEP=(DTHE(N)-DTHE( 1) )/K2 

WRITE(6,302) 
302    FORMAT(«  •  ,4X, ' THET ' ,  12X, ■ EMAG *  ) 

DC  3  K=1,K2 

THETA(K)=THET 

EN  =  A(  1) 

DC  4  J=2,L 
4      EM=EM+A( J)*THETA(K)**< J-l) 

EPAG(K)=EM 

WPITE(6,301)  THET,EMAG(K) 

301  FORMATP  '  ,  F  10.  5  ,  5X  ,F  1  0  .7  ) 
TFET=THET+STEP 

3      CONTINUE 

C      CALCULATE  THE  APERATURE  FIELD  STRENGTH 

DC  7  K=1,K2 

ENA(K)  =EMAG(K)*SQRT( 1 . 0/COS ( THET A ( K ) ) ) 

WR1TE(6,301)THETA(K) ,EMA(K) 

7  CONTINUE 

DO  8  K=1  K? 

TFETA(K )=THETA(K)*180.C/3.1415 

8  CONTINUE 

CALL  DRAW(K2  , THET A , EMAG , I , 0 , L A E EL ,T I TL E , 2C .0 ,0 , 0 , 3 , 2, 
12,6,7, I, LAST) 

CALL  DRAW(K2  ,THETA,EMA  ,3  ,0  ,  LAE  EL  ,T  ITL  E  ,  2G  .CO,  0  ,  3,  2, 
12,6,7, 1,LAST) 

STCP 
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END 
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PHE-CCMPCNENT  LUNEBERG  LENS  PATTERNS 


C 
C 

c 
c 
c 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 

c 

c 
c 
c 
c 
c 

12 

4 


5 

c 
c 

100 


c 
c 

13 


THIS  PR 
OF  THE 
LENS  WI 
PCLAR1Z 
PCLARIZ 
IMPLICI 
COMPLEX 
M  =  C 

STEP  IS 
THE  FIE 
THETAI 
FIELD  S 
THETAF 
FIELD  S 
ANGPD  I 
IN  WHIG 
DEGREES 
O-DEGRE 
EC. 

READ! 5, 
FCRMAT( 
PIE=3.1 
ANGP=AN 
CALCULA 
NUMBER 
NUM=(TH 
MC  IS  T 
TC  BE  P 
RADIUS 
FREQ  IS 
CALCULA 
THICK  I 
READ(5. 
FORMAT( 
IF (MO  . 
N  =  N+1 
WRITEtc 
1THICK,A 
FORMAT ( 
1  «,'STE 
2»FREQUE 
3, 'THICK 
R=200.0 
AK  IS  T 
AK=20.0 
CJKR=(C 
CJKRN  I 
CJKRN=C 
WRITE(6 
FORMAT ( 
CALL  EP 
1,NC  ,FR 
EPHEE  W 
OF  THE 
I F ( MO  . 
STCP 
END 


OGRAM  WAS  DESIGNED  TO  CALCULATE  PF E-COMFONENT 
RADIATION  PATTERNS  OF  THE  CYLINDRICAL  LUNEBERG 
TH  3EAM  CENTERED  AT  O-DEGR E ES , FOR  A  VERTICALLY 
ED  FEED. THIS  PATTERN  IS  ICENTICAL  FCP  ANY  FEEC 
ATION  IN  THE  X-Y  PLANE 
T  COMPLEX (C) 
CJ/(0.O, 1.0)/ 

THE  STEP-SIZE  IN  DEGREES  BETWEEN  PCINTS  WHERE 
LD  STRENGTH  IS  CALCULATED. 

IS  THE  INITIAL  ANGLE  IN  DEGREES  FCP  WHICH  THE 
TRENGTH  IS  CALCULATED. 

IS  THE  FINAL  ANGLE  IN  DEGREES  FOR  WHICH  THE 
TRENGTH  IS  CALCULATED. 

S  THE  ANGLE  IN  CEGREES  THAT  CONTROLS  THE  PLANE 
H  THE  LENS  PATTERNS  ARE  CALCULATED,  IF  90 

THE  THE  E-PLANE  PATTERNS  ARE  CALCLLATED,  IF 
ES  THEN  THE  VERTICAL  BEAM  PATTERN  IS  CALCULAT- 


2)  ST 

F12.3 
41592 
GPD^{ 
TE  NU 
OF  PO 
ETAF- 
HE  IN 
LOTTE 
IS  Tri 

THE 
TED. 
S  THE 
4)  NO 
Il/Fl 
EQ.O) 


EP,T 
/F12 

fc53!> 
PIE/ 
MBER 
INTS 
THET 
TEGE 
D  BY 
E  RA 
FREG 


HETAI,  THETAF, ANGPC 
.3/F12.3/F12.3) 


18C.0) 

OF  STEPS  ,NUM, 

TC  BE  PLOTTED 
AI  J/STEP  +1 
R  THAT  CONTROLS  THE  NUMBER 

DRAW   IN  SLBROLTINE  EPhEE. 
D1US  OF  THE  LUNEBERG  LENS. 
UENCY  FOR  WHICH  THE  PATTERNS 


WHICH  IS  ALSO  THE 
BY  SUBROUTINE  CRAW 


CF  PLOTS 


ARE 


HEIGHT  OR  X-DIMENSICN 
,  RADIUS,  FREQ, THICK 

2.3/F12.3/F12.3) 

GO  TO  14 


OF  THE  LENS 


,5)THETAI,  THETAF, STEP,  R AD  1  US , FR EC  ,  NLM , MO  , 

NGPC 

•  '  j'THETAI^jFll.b,/,'  •  , 'ThETAF='  ,F11.6,/,  • 

P=»,FU.6,/,«  •  ,  'RADIUS='  ,F11.6,/ ,  •  ', 

NCY  =  «  ,FU  .6,/,  •  •  ,  «NUM=»  ,  16,/,  •  ','NC   ',16,/' 


=  •  ,FU  .6,/,  • 


•PHE=» ,F I  1.6) 


HE  PRO 

*PIE*F 

J*AK^R 

S  THE 

JKR/CA 

,100)C 

'  •  ,'C 

HEE(AK 

EG,M) 

ILL  CA 

E-F1EL 

EQ.CJG 


PAGATICN  CCNSTANT. 

REQ/3.0 

ADIUS/(4.0*PIE*R) )*CEXP(CJ*AK*(RADIUS-R)) 

NORMALIZED  COEFFICIENTS 

BS(CJKR) 

JKRN 

JKRN=« ,F  1  l.9,2X,» J* ,F11.9) 

, THICK, AN GP, THETAI, STEP,  NUM,CJKRN,RACILS 

LCULATE  AND  PLOT  THE  PHE-CONPCNENTS 
D  OF  THE  LUNEBURG  LENS 
0  TO  13 


SLBROUTINE  EPHEE  ( AK  1 , THI , ANGP 1 , ANGT I , STE , NLM ,C JKRN , 
'IRA,  MOD,  FREi  Ml) 
C      THIS  SUBROUTINE  WAS  WRITTEN  TO  CALCULATE  ANC  PLOT 
C      THE  PHE-COMPONENT  OF  THE  ELECTRIC  FIELD  CF  THE  LUNE- 
C      BERG  LENS  FOR  ANY  FEED  FUNCTION  CR  ANY  MEASURED  FEED 
C      DISTRIBUTION  THAT  CAN  BE  APPROXIMATED  BY  A  CURVE 
C      FITTING  TECHNIQUE. 

IMPLICIT  COMPLEX (C) 
C      THE  NEXT  FOUR  STATEMENTS  CONTROL  THE  LABELING  AND 
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LOO 


C 
C 
C 

c 
c 
c 
c 


c 
c 
c 


900 

1000 
2 


I 

c 
c 


3 

84 


85 


23 
C 


TITLING 

REAL  LA 
REAL  *3 
READ( 5, 
FORMAT ( 
COMPLEX 
EXTERNA 
CFCT1  I 
LUNEBER 
DISTRIB 
CFCT3  I 
LUNEBER 
D  1STRI3 
ANGULAR 
DIMENSI 
1FREQ15) 
CCMMON/ 
COMMON 
ThIOTH 
RAD=RA 
RADHM1 
FREQ(Ml 
AK=AK1 
ANGP=AN 
PIE=3.1 
THET=AN 
TH2=-TH 
TH3=THI 
ANCRM=0 
A NORM  I 
FIELD  S 
ENTER  L 
DO  1  1  = 
L1=I 
L3  =  I 
THETR=T 
CALL  CW 
Ll=2 
XL=ALIM 
XU=ALIM 
CALL  CW 
IF(L3.N 
L3=2 
GO  TO  9 
CEPH=CE 
IF(ANGR 
AN0PM=C 
THETA(  I 
CEPHE(  I 
THET=TH 
CCNTINU 
ENTER  L 
OF  THE 
DC  3  J  = 
CEPHE( J 
EPHEM  (J 
CCNTINU 
WRITE16 
FORMAT 
1' RADIUS 
2,  'THETA 
3'STRENG 
J  =  0 

DO  23  I 
WRITE  ( 
FORMAT! 
J=J  +  l 
IF(J.LT 
WRITE ( 6 

J=0 

CCNTINU 

ENTER  L 


OF  THE  G 
BEL  /' 

riTLt(L2 
100)  TITL 
6A3) 

CJ/(0.0, 
L  CFCT1,C 
S  A  FUNCT 
G  LENS  IN 
UTION  IN 
S  A  FUNCT 
G  L F N S  IN 
UTION    AS 

DISPLACE 
ON  THETA( 

ALIM/ALIM 
AK»  ANGP,T 
I 


RAPHS  PLOTTED  BY  SUBROUTINE  CRAW. 

'/ 
) 
E 


I  .0) 
FCT3 
ION 

TER 
THE 
ION 

TER 
A  FU 
ME  NT 
40C) 


/ 


THE 


DESCRIBING  THE  FAR  FIELD  OF 

MS  OF  THE  APPERATURE  FIELD 

X-DIRECTION. 

DESCRIBING  TFE  FAR  FIELC  OF 

MS  OF  THE  APPERATURE  FIELD 

NCTION  OF  THE 

FROM  BEAM  CENTER. 
,CEPHE(400)  ,EPHEM(400),RADI(5)t 


THE 


1 ,ALIM2,THIC 
HETR»RADtL3 ,L1 


)=RA 
)  =  FRE 

GP1 

415926535 

GTI 

1/2.0 

/2.0 

.0 

S  A   VARI 

TPENGTH. 

OOP  TO  CA 

l,NUM 


FET*(PIE/ 
EDF(CFCTl 

L 
2 

EDF(CFCT3 
E.-l)  GO 

00 

PHl*CEPH3 

M.GE.CABS 

ABS(CEPH) 
)  =  THET 
)=CEPH 
ET+STE 

COP  TO  NO 
E-PIELD  ( 
1,NUM 
)=CEPHE( J 
)=CA6S(CE 


ABLE    USED    TO    NORMALIZE    THE    ELECTRIC 
LCULATE    THE    E-FIELO    FOR    NUM-POINTS 


130.0) 

,TH2,TH3, 15.CEPH1) 


,     XL     ,    XU     ,375,CEPh3) 
TO     1000 


(CEPH) )    GO    TO    2 


,  'PHE 
t  F  3  . 4 


,34)RADI  ( 
<  •  1« 
OF' 
' ,  12X,  »FI 
TH'  ) 

=1, NUM 
6,85)THET 
'     '  ,F11.3 

.72  )G0    TO 

,S4)RAD,P 


RMALIZE    AND    CALCULATE    THE    MAGNITUDE 
EPHEM  )  . 

) /ANORM 
PHE(J)  ) 

Ml) , FREQ(Ml) 

-COMPONENT  OF  THE  E-FIELD  FCR  LENS', 
,2X, 'METERS  AT • , F8 . 3 , 2X , « GFZ ■ , // ■  '  ,  6X 
ELD  STRENGTH*  ,10X, 'COMPLEX  FIELD  ', 


A(  I  )  ,EPHEM(  I  ),CEPHE(  I) 
,8X,F14.6,12X,F10.6,F10.6) 

23 

RE 


COP  TO  CALCULATE  THE  FIELC  MAGNITUCES  IN  DECI- 
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C  BELS. 

DC    20    1=1, NUM 

EPHEM( I )=20. 0*AL0G10( 5PHEM( I ) ) 
20  CCNTINUE 

CALL    DRAW(NUM,THETA, EPHEM , MOD, 0 , LA  EEL , T  IT L E , 0, 
12,2,6,3, 1,LAST) 
90  RETURN 

END 


0,7,2, 


C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 


500 


700 
800 


900 
901 

1000 


SUBROUTINE  CWEOF  ( CF, XL , XU, NX, C ANSI 

CWEDF  WAS  OBTAINED  FROM  LCD4.  E.G.NEELY  . 

CWEOF  IS  A  SUBROUTINE  WHICH  WILL  NUMERICALLY 
INTEGRATE  A  USER  SUPPLIED  FUNCTION  BETWEEN  SPECIFIED 
LIMITS. (SINGLE  PRECISION) 

CF    -  NAME  OF  FUNCTION  SUBPROGRAM.   NUST  BE 
LISTED  IN  AN  EXTERNAL  STATEMENT. 

XL    -  LOWER  LIMIT  OF  INTEGRATION 

XU    -  UPPER  LIM!T  OF  INTEGRATION 
NX    -APPROXMATE  NUMBER  OF  NODES  AT  WHICH  TO 
EVALUATE  THE  FUNCTION. 

CANS  -  RESULT  OF  INTEGRATION 

IMPLJ.CI  T  COMPLEX"^  (C) 
REAL*8  DXDX,XX 

REAL*4  CW(6)/3  2.,216.,27.,272.,27.,216./ 
CCMMON  AK,ANGP,THETR,RAD,L3 ,L1 
IF(NX.LE.O)  GO  TO  900 
N=( (NX+4)/6)*6+l 
DX=(XU-XL)/FL0AT(N-1) 
DXCX=DBLE(DXJ 
NWIX=N/6 
X=  XL 

CANS=-CF(X)*4i.O 
IFCL3.EC.-1J  GO  TO  1000 
IF(Ll.EO.l)  GO  TO  500 
IFU3.NE.1  J  GO  TO  500 
L2  =  2 

DO  8C0  MX=l,NWIX 
DC  700  KX=1,6 
CANS=CANS+CW(KX)*CF(X) 
XX=DBLE (X) 
X=SNGL(XX+DXDX) 
CONTINUE 
CONTINUE 

CANS=(CANS+41 .0*CF(X)  )*DX/ 140.0 
RETURN 

WRITE (6 ,901 J  N 
FCRMAT( • OERROR  IN  CALLING  PARAMETER  *****   N  =  », 


RETURN 
END 


ii^^ijt^  •  •//  ) 


c 
c 
c 


FUNCTI 

CFCT1 

VARIAT 

BERG  L 

IMPLIC 

CCMPLE 

CCMMON 

COMMON 

CEX=CJ 

CFCT1= 

RETURN 

END 


DESCRIBES  THE 
AS  A  FUNCTION 


ON    CFCTKX  ) 

IS    A    FUNCTION    WHICH 

IONS    OF    THE     E-FIELD 

ENS    HEIGHT. 

IT    COMPLEX (C) 

X    CJ/( 0.0,  1.01/ 

AK,ANGP,THETR,RAD,L3,L1 
/ALIM/ALIMltALI M2,THIC 
*AK*SIN(THETR)*COS( ANGP)*X 
CEXP(CEX)*C0S(3. 14159* X/THIC) 


FAR-FIELD 
OF    THE    LUNE- 
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c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


200 
400 
60 
50 


6 
300 

600 

4 
C 
C 


C 
C 

c 
c 
c 
c 
c 


FUNCTTO 
CFCT3  I 
VAX  I  AT  I 
(THE  AN 
UTILIZE 
AND  COM 
OF  THE 
IMPLICI 
INPLICI 
CCMPLEX 
D I  MENS  I 

LDELY(30 
REAL*3 
CCMMON 
CCMMON/ 
L3  WILL 
AND  WIL 
IN  ORDE 
DESCRIB 
IF  L3  I 
CALCULA 
IF(L3.N 
READ(5, 
FCPMAT( 
READ( 5, 
FORMAT ( 
DC  60  J 
DMAG( J  ) 
DO  50  I 
DTHE( I) 
ALIML=D 
ALIM2=0 
DO  5  1  = 
DW I (  I  )  = 
CALL  LS 
L=1+IAB 
DC  6  1  = 
A(I  )=DA 
THETA( I 
WRITE(6 
FORMAT { 
L3=-l 
RETURN 
E  f  =  A  ( 1  > 
DO  4  J  = 
EM=EM+A 
EA1  IS 
OVER  TH 
EA1=EM 
CER=CJ* 

1SI.N(THE 
THE  SQR 
EQUATIO 
THE  REM 
FIELD  M 
IFY  THI 
FIELD  P 
TERMS  A 
CFCT3=C 

LSINCANG 
RETURN 
END 


N    CFCT3(ANG 
S    A    FUNCTIO 
ONS    OF     THE 
GLE    FROM     BE 
D    TO    FIT    A 
BINES    THIS 
FAR-FIELD    E 
T    REAL*8(D) 
T    COMPLEX     ( 

CJ/(0.0,  1  . 
ON  DTHETAd 
), A(20) ,DSA 
DTITL(  10J/1 
AK,ANGP,THE 
ALIM/ALIM1, 

BE    EQUAL     T 
L    PROCEED    T 
R    TO    CALCUL 
E    THE    FEED 
S    NOT     L,     TH 
TED    BY    A    PR 
E. L)    GO    TO 
20  0)    M,N 
13/  12) 

400)  (DTHET 
F10.7,F12.6 
=  1,N 

=DSQRT(DMAG 
=  1,N 

=DTHETA(I )* 
THE( I) 
ThE(N) 
UN 
L.ODO 

CPL2(N,M,DT 
S(M) 
ltL 
(I) 

)=DTHETA(I ) 
,300)  (I,A( 
*     •  ,'  AC  ,  14 


2,L 

(J )*ANG**( J 
A    FUNCTION 
E    PLANE    OF 


) 

N  WHICH  DESCRIBES  THE  FAR 
E-FIELD  AS  A  FUNCTION  OF 
AM  CENTER).  TO  DO  THIS  LS 
CURVE  TO  THE  MEASURED  FEE 
WITH  THE  OTHER  ANGULAR  VA 
QUATIONS  . 

C) 
0)/ 

00)  ,DMAG(  100)  tDA(20) , DW  I  ( 
(  20)  ,THETA(  10  J)  , EMAG(  LOO 
0*'  •/ 
TR,RAD,L3,L1 

ALIM2,THIC 

C  1  ON  THE  FIRST  CALL  OF 

HPOUGH  THE  NEXT  STEPS  DCW 

ATE  THE  POLYNOMINIAL  COEF 

FUNCTION. 

EN  THE  POLYNOMINIAL  COEF. 

EVIOUS  PASS. 

600 


A(I  )  ,DMAG(I)  ,1=1, N) 
) 

(J)  ) 

(3. 14L59265D0/180.D0) 


HE   ,  DMAG,DWI,CY,DELY,DA,DSA,DT  ITL) 


FIELD 
ANG 

CPL2  IS 
D  FUNCTION 
PI ATIONS 


30),DY(80), 
,DTHE( LOO) 


CFCT3 

N  TO  600 

.  THAT 

HAVE  BEEN 


I  ),  1=1, L) 

»'  )  =  •  ,FL8.6) 


-L) 

WHICH  DESCRIBES 
THE  LENS. 


THE  FEEC  VARIATIONS 


AK* 
TR) 
T(l 

N  A 

AIN 

ETH 

S  P 

ATT 

RE 

EXP 

P)* 


RAD*(COS 
*SIN( ANG 
,0/COS (A 
SSOCI ATE 
OER  IS  F 
GD  TO  TH 
ROGRAM  T 
ERNS  CHA 
ALREADY 
(CER)*<S 
(COS(THE 


(ANG 
P)  ) 
NG)  ) 

D  WI 
ROM 
E  SE 
0  PR 
NGE 
INCL 
IN(T 
TR)  + 


)*(C0S(THETR)-1.0)-SIMANG)* 

T  OF  THIS 


TER 
TH  T 
THE 
MI-C 
EDIC 
THIS 
UDED 
HETR 
1.0) 


M  IS  THE  GNLY  PAR 
HE  LENS  FUNCTICN. 
APPL ICATION  OF  TH 
IRCULAR  APERATURE 
T  THE  CCNSTANT-K 
FUNCTION.  THE  CO 

)*SI.N(  ANG)-COS(  AN 
)*SQRT(1  .0/CQ5( AN 


E  APERATURE 
.TO  MOD- 
LENS  FAR 
RRECT  PHASE 

G)* 

G)  )*EAL 
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THETA-COMPONENT  LUNEBERG  LENS  PATTERNS 

C      THIS  PROGRAM  WAS  DESIGNED  TO  CALCULATE  TEE  THETA  COM- 
C      PONENT  OF  THE  CYLINDRICAL  LUNEBERG  LENS  RACIATION  PAT- 
C      TERNS  WITH  BEAM  CENTERED  AT  O-D EGR EES , FQ R  ANY  FEED 
C      POLARIZED  IN  THE  X-Y  PLANE. 

IMPLICIT  COMPLEX(C) 

COMPLEX  CJ/(0.0,1 .0)/ 

M  =  0 
C      STEP  IS  THE  STEP-SIZE  IN  DEGREES  BETWEEN  PCINTS  WHERE 
C      THE  FIELD  STRENGTH  IS  CALCULATED. 

C      THETAI  IS  THE  INITIAL  ANGLE  IN  DEGREES  FOR  WHICH  THE 
C      FIELD  STRENGTH  IS  CALCULATED. 

C      THETAF  IS  THE  FINAL  ANGLE  IN  DEGREES  FOR  WFICF  THE 
C      FIELD  STRENGTH  IS  CALCULATED. 

C      ANGPD  IS  THE  ANGLE  IN  DEGREES  THAT  CONTROLS  THE  PLANE 
C      IN  WHICH  THE  LENS  PATTERNS  ARE  CALCULATED,  IF  90 
C      DEGREES  THE  E-PLANE  PATTERNS  ARE  CALCULATED  ,  I F 
C      O-DEGREES  THEN  THE  VERTICAL  BEAM  PATTERN  IS  CALCULAT- 
C      ED. 

READ(5,2>  STEP,THETAI,  THETAF, ANGPD 
2      F0RMAT(F12.3/F12.3 ,F12.3/F12.3) 

PIE=3.L415926535897 

ANGP=ANGPD*( PIE/ 18  0.0) 
C      CALCULATE  NUMBER  OF  STEPS  REQUIREC  WHICH  IS  THE  SAME 
C      AS  THE  NUMBER  OP  POINTS  TO  BE  PLOTTED  BY  DRAW 

NUM=(THETAF-THETAI )/STEP  +L 
C      MO  IS  THE  INTEGER  THAT  CONTROLS  THE  NUMBER  CF  PLOTS 
C      TO  BE  PLOTTED  BY  DRAW   IN  SUBROUTINE  ETHETA 
C      RADIUS  IS  THE  RADIUS  OF  THE  LUNEBERG  LENS. 
C      FREQ  IS  THE  FREQUENCY  FOR  WHICH  THE  PATTERNS  ARE 
C      CALCULATED. 

C      THICK  IS  THE  HEIGHT. OF  THE  LENS. 
L2     READ(5,4)  MO  ,  RADIUS,  FREG, THICK 

4  FCRMAT(  I1/F12.3/F12.3/F12.3) 
IF(MO  .EQ.0)GO  TO  14 

M=M+1 
14     VvPITE(6,5)THETAI  ,  THETAF, STEP,  RADIUS  ,  FREQ,  MJM,  MO, 
1THICK, ANGPD 

5  FORMAT (■  ' , •THETAI=,,F11.3, /, '  » ,  'THETAF=  '  , F  1 1. 3, / , 
1«  «  ,' STEP='  ,F11.6, /,•  •  ,' RADIUS  =  '  ,F11.6,/,  '  », 
2»FREQUENCY  =  « ,F11.6,/,  •  • ,  • NUM= •  ,  I  6  ,  / , ■  'j'NC  =,',16 
3,/'  •  ,'THICK=« ,F11  .6,/'  ' ,  'PFE=',F11.6) 

R=200.0 
C      AK  IS  THE  PROPAGATION  CONSTANT. 

AK=20.0':PlE*FREQ/3  .0 

CJKR=(CJ*AK*RADIUS/(4.0*PIE*R) ) *CEXP (C J* AK* ( R AD I US-R ) ) 

CJKRN=CJKR/CABS  (CJKR) 

WRITE( 6,100)CJKPN 
100  FORMATt1     •  ,«CJKRN=  •  ,F  ll.9,2X, ' J '  ,F11.9) 

CALL    ETHETA( AK, THICK, ANGP, THETA I, STEP, NUM , C JKRN ,R ADI US 
IMC     ,FREQ,M) 
C  EFTETA    WILL    CALCULATE     AND    PLOT    THE    THETA    COMPONENTS 

C  -FIELD    PATTERNS    OF    THE       LUNEBURG    LENS 

IF(MO    .EQ.0)GO    TO    13 
13  STCP 

END 

SUBROUT INE  ETHETA ( AK 1 , THI , ANGP 1 , ANGTI ,STE, NUM, C JKRN, 
1RA,M0D,FRE,M1 ) 
C      THIS  SUBROUTINE  WAS  WRITTEN  TO  CALCULATE  ANC  PLOT 
C      TFE  THETA-COMPONENT  OF  THE  ELECTRIC  FIELD  OF  THE  LUNE- 
C      BERG  LENS  FOR  ANY  FEED  FUNCTION  OR  ANY  MEASURED  FEED 
C       DISTRIBUTION  THAT  CAN  BE  APPRI X I PATED  BY  A  CURVE 
C       FITTING  TECHNIQUE. 

IMPLICIT  CCMPLEX(C) 
C      THE  NEXT  FOUR  STATEMENTS  CONTROL  THE  LABELING  AND 
C      TITLING  OF  THE  GRAPHS  PLOTTED  BY  SUBROUTINE  DRAW. 

REAL  LABEL  /•      • / 
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100 

c 
c 
c 

c 
c 

c 
c 


900 


1000 


96 

c 

15 


95 

1 

C 

c 


3 
84 


85 


REAL  *8 
PEAD(5, 
FORMAT! 
COMPLEX 
CFCT1  I 
LUNEBER 
DISTRIB 
CFCT4  I 
LUNESER 
DISTRIB 
ANGULAR 
EXTERNA 
DIMEKSI 
1FREQ(5) 
COMMON 
COMMON/ 
THIC=TH 
R£D=RA 
R  A  C I  ( M 1 
FREQ(M1 
AK  =  AK1 
ANGP=AN 
PIE=3.1 
TFET=AN 
TF2=-TH 
TH3=THI 
ANORM=0 
ENTER  L 
DO  1  1  = 
L  1=1 
L3  =  I 
THETR=T 
CALL  CW 
Ll=2 
XL=ALIM 
XU=ALIM 
CALL  CW 
IF(L3.N 
L3  =  2 
GC  TO  9 
CETH=CE 
WRITE (6 
FORMAT! 
WRITE(6 
FORMAT! 
IF(ANOR 
ANCRM=C 
THETAd 
CETHE! I 
TFET  =  TH 
WRITE!6 
FORMAT! 
CCNTINU 
ENTER  L 
OF  THE 
DC  3  J  = 
CETHE! J 
ETHEM(J 
CCNTINU 
WRITE! 6 
FORMAT 
1' RADIUS 
2,  'THETA 
3' STRENG 
J  =  0 

DO  23  I 
WRITE  ( 
FORMAT! 
J  =  J  +  1 
IFtJ.LT 
WRITE ( 6 
J=0 


I  ELD    OF 
FIELD 

IELD    OF 
FIELD 


TITLE (121 

10  0)    TITLE 

6A8) 

C  J / ( 0 . 0,1.0)/ 
S  A  FUNCTION  DESCRIBING  THE  FAR  p 
G  LENS  IN  TERMS  OF  THE  APPERATURE 
UTION  IN  THE  X-DIRECTION. 
S  A  FUNCTION  DESCRIBING  THE  FA}  F 
G  LENS  IN  TERMS  OF  THE  APPERATURE 
UTION    AS     A    FUNCTION    OF    THE 

DISPLACEMENT    FROM    BEAM    CENTER. 
L    CFCT1TCFCT4 
ON    THETA (400) , CETHE (400) , ETHEM! 400  )  , RAD I ( 

AK, ANGP,THETR,RAD,L3,Ll 
ALIM/ALIM1 , ALIM2,THIC 
I 

)  =  RA 
)=FRE 


GP1 

41592653589 

GTI 

1/2  .0 

/2.  C 

.0 

COP    TO    CALCULATE 

1,NUM 


THE 


THE 


5) 


THE    E-FIELD    FOR    NUM-POINTSIN 


HET*( PIE/180.0) 
EDF(CFCTL,TH2,TH3t 5»CETH1) 


1 

2 

EDF(CFCT4,     XL        , 

E.-l)    GO    TO    1000 


XU        ,375,CETH4) 


00 

TH1 

,96 
<     i 

,  15 
■     i 

M.G 
A8S 
)=T 
)=C 

ET  + 

,95 
i     i 


*CETH4 

)    CETH1,CETH4,CETH 

,6E12.6) 

JCETH 

,'CETH=«  ,2F10.5) 

E.CABS(CETH) )    GO    TO    2 

(CETH) 

HET 

ETH 

STE 

)     TF.ET,ANCRM 

,2F12.5) 


OOP 
E-FI 
1,  NU 
)=CE 
)  =  CA 
E 

,34) 
('  1« 
OF' 
1  tl2 
TH»  ) 


AND    CALCULATE    THE    N4GNITUCE 


TO    NORMALIZE 

ELD    (EPHEM). 

M 

THE( J) /ANORM 

BS(CETHE( J)  ) 

RADI ( Ml)  ,FREQ(M1 > 

, 'THETA-CGMPONENT    OF     THE    E-FIELD    FOR    LENS' 

, F8. 4, 2X, 'METERS    AT  '  ,  Fl 

X,  'FI ELD  STRENGTH'  ,10X 


t  <1A  , 

•COMPLEX 


FIELD  • 


=1, NUM 

6,85)THETA(  I )  ,  ETHEM (  I),CETHE(  I) 

'  «  ,F11.3,8X,F14.6,12X,F10.6,F10.6) 

.72  )GC  TO  23 
f84)RA0TFRE 
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23     CONTINUE 

C      ENTER  LOOP  TO  CALCULATE  THE  FIELC  MAGNITUDES  IN  DECI- 

C      BELS. 

DC  20  I =1 ,NUM 

ETHEM( I )=2  0. 0*AL0G10(ETHEM(  I)  ) 
20     CONTINUE 

CALL  DRAW(NUMtTHETA, ETHEM,MOD,0,LAEEL,TITLE,0,0,7,3,2, 
12,6,8, L,LAST) 
90     RETURN 

END 

SUBROUTINE  CWEDF  ( CF , XL , XU, NX, C ANS ) 
C 

C      C^EDF  WAS  OBTAINED  FROM  LCDR.  E.G.NEELY  . 
C         CWEDF  IS  A  SUBROUTINE  WHICH  WILL  NUMERICALLY 
C      INTEGRATE  A  USER  SUPPLIED  FUNCTION  BETWEEN  SPECIFIED 
C      LIMITS. (SINGLE  PRECISION) 

C  CF    -  NAME  OF  FUNCTION  SUBPROGRAM.   MUST  BE 

C      LISTED  IN  AN  EXTERNAL  STATEMENT. 
C  XL    -  LOWER  LIMIT  OF  INTEGRATION 

C  XU      UPPER  LIMIT  OF  INTEGRATION 

C  NX    -APPROXMATE  NUMBER  OF  NCDES  AT  WHICH  TO 

C      EVALUATE  THE  FUNCTION. 
C  CANS  -  RESULT  OF  INTEGRATION 

C 

IMPLICIT  COMPLEXES  (C) 

REAL* 3  DXDX,XX 

REAL*4  CW( 6) /8  2. ,216. ,27.  , 272 . , 27  .  ,  216  .  / 

CCMMON  AK, ANGP,THETR,RAD,L3,L1 

IF(NX.LE.O)  GO  TO  900 

N=( (NX+4)/6)*6+l 

DX=(XU-XL)/FLOAT (N-l) 

DXDX=DBLE(DX) 

NWIX=N/6 

X  =  X  L 

CANS=-CF(X)*4l.O 

IF(L3.EQ.~1)  GO  TO  1000 

IF(L1 .EC.l)  GO  TO  500 

1F(L3.NE. 1)G0  TO  5C0 

L3=2 
500    DC  800  MX=1 ,NWIX 

DO  700  KX=1,6 

CANS=CANS+CW(KX)*CF(X) 

XX=DBLE(X) 

X=SNGL( XX+DXDX) 
700   CONTINUE 
800   CONTINUE 

CANS=(CANS+41.0*CF(X) )*DX/140.0 

RETURN 

900  WRITt(6,90i)  N 

901  FORMAT ( 'OERROR  IN  CALLING  PARAMETER  *****   N  =  ', 
L I S  , «   ****** // ) 

1000   RETURN 

END 

COMPLEX  FUNCTION  CFCTKX) 
C      CFCT1  IS  A  FUNCTION  WHICH  DESCRIBES  THE  FAR-FIELD 
C      VARIATIONS  OF  THE  E-FIELD  AS  A  FUNCTION  OF  THE  LUNE- 
C      BERG  LENS  HEIGHT. 

IMPLICIT  COMPLEX(C) 

COMPLEX  CJ/CO.O, 1.0)/ 

COMMON  AK,ANGP,THETR,RAD,L3  ,Ll 

COMMON/ AL I M/ AL I M I , AL I M2 , TH I C 

CEX=CJ*AK*SIN(THETR)*COS(ANGP)*X 

CFCT1=CEXP(CEX)*C0S(X*3.14159/THIC)*(1 . O+CCS (THET R ) )* 
lCCS(ANGP) 

RETURN 

END 

COMPLEX  FUNCTION  CFCTA(ANG) 
C      CFCT4  IS  A  FUNCTION  WHICH  DESCRIBES  THE  FAR  FIELD 
C      VARIATIONS  OF  THE  E-FICLD  AS  A  FUNCTION  CF  ANG 
C      (THE  AMGLE  FROM  BEAM  CENTER).  TO  DO    THIS  LSuhL?.  IS 
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